INDEPENDENT SETS IN CERTAIN CLASSES OF (ALMOST) 

REGULAR GRAPHS 



Alexander Burstein 
Department of Mathematics 

Iowa State University 
Ames, IA 50011-2064, USA 
burstein@math . iastate . edu 

Sergey Kitaev 
Department of Mathematics 
University of Kentucky 
Lexington, KY 40506-0027, USA 
kitaev@ms . uky . edu 

Toufik Mansour 
Department of Mathematics 
Haifa University 
31905 Haifa, Israel 
toufik@math.haif a. ac . il 

Abstract. We enumerate the independent sets of several classes of regular and al- 
most regular graphs and compute the corresponding generating functions. We also note 
the relations between these graphs and other combinatorial objects and, in some cases, 
construct the corresponding bijections. 

Keywords: independent sets, regular graphs, transfer matrix method 

2000 Mathematics Subject Classification: Primary 05A05, 05A15; Secondary 
30B70, 42C05 

1. Introduction 

Let 1(G) denote the number of independent sets of a graph G. This number can be 
determined for some special classes of graphs (see [7j for a survey). For instance, 1(G) 
was studied for grid graphs (see pQ), multipartite complete graphs, and path- and cyclic- 
schemes. Many of these numbers are given by certain combinations of Fibonacci numbers, 
some others by Lucas numbers. 

In this paper, we study four classes of graphs. To define these classes, we recall that 
a line graph L(G) of a graph G is obtained by associating a vertex with each edge G 
and connecting two vertices with an edge if and only if the corresponding edges of G are 
adjacent. Also, recall that a cycle graph Ce is a graph on i nodes containing a single cycle 
through all nodes. 
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We now give the definitions of our classes. 

Class 1: Let G\ = Ci, the £-cycle graph. We obtain Gj by superimposing the 
line graph L(G\) onto the graph G\, that is splitting each edge of Gj with the 
corresponding vertex of L{G\) and then adding the edges of L{G\). More generally, 
Gg is obtained by superimposing L n ~ l (G\) onto G™~ . For example, in Figure [TJ 
we have the graphs G\ and G\, respectively, if one ignores the dashed edges. 
Clearly, all but n of the £ n nodes of G™ have degree 4, and we say that G™ is an 
almost 4-regular graph. 




Figure 1. Examples of (almost) 4-regular graphs under consideration. 



Class 2: The graph R™ is obtained from G™ by duplicating the edges of G\. For 
example, in Figure ^ the extra edges are the dashed edges, and by adding them 
we get the graphs R\ and R\ respectively. So, to get R™ we add £ additional edges 
to G"', and it is easy to see that R% is a 4-regular graph. 

Class 3: Let K\ = Kg, a complete graph on i nodes. Put the £ nodes of K\ on a 
circle and draw the remaining (£ — 1)! — £ edges. Call the first I edges external 
and the remaining edges, internal. Then construction of K™ is similar to that of 
G™, except that: 

(1) The basis of the construction is now Kg, rather than C#. 

(2) On each iteration i, the graph superimposed onto R\ _1 is not L % {K() but 
rather the complete graph on the nodes of L l (Ci), the line graph of Ci formed 
by the external edges of K\. 

In Figure El we show how to construct K\ from K\ (the dashed edges should 
be ignored). In that figure, the external edges of respective complete graphs 
are in bold. We also remark that the internal edges do not intersect each other. 
Moreover, Figure El suggests a convenient way of representing Kf, where each node 
of the graph lies only on internal edges incident with that node. We achieve that 
by putting the nodes of the line graphs under consideration off the centers of the 
corresponding edges. Indeed, if we were using the centers of the external edges, 
K\ would look as in Figure El which is misleading since, for example, the node a 
does not lie on the edge be. 

The graph Kf is almost (£+l)-regular, since all but £ of its £ n nodes have degree 
£ + 1. Moreover, it follows from our definitions that = G%. 




FIGURE 2. A way of constructing K\ and P 4 3 . 

Class 4: The graph P" is obtained from the graph by duplicating the external 
edges in the graph K\. For example, in Figure El the extra edges are the dashed 
edges, and by adding them we construct the graph P| from K\. So to get P™ we 
add i extra edges to Kg, and it is easy to see that Pf is an (£ + l)-regular graph. 



b 




c 



Figure 3. A bad presentation of K\. 



Let gt{n), Ti(ri), kf(n) and pe(n) denote the number of independent sets in the graphs 
G", R™, K™ and P™ respectively. In our paper we study these numbers. In Section EJ we 
give an algorithm for calculating all these numbers. For the numbers pe{n), we provide 
an explicit generating function (see Theorem I3.4|) . However, in order to illustrate our 
approach to the problem, in Section |21 we consider g^{n) and find an explicit formula for 
it. 

Our choice of the graphs to study was motivated by the so called de Bruijn graphs, 
which are defined as follows. A de Bruijn graph is a directed graph G n = G n (V, E), where 
the set of vertices V is the set of all the words of length n in a finite alphabet A, and 
there is an arc from Vi = (vn, . . . , v in ) to Vj = (vjx, . . . , vj n ) if 

Vi2 = Vji, V i3 = V j2 , ...,V in = Uj(n-1), 

that is when the words Vi and Vj overlap by (n — 1) letters. 

The de Bruijn graphs were first introduced (for the alphabet A = {0, 1}) by de Bruijn 
in 1944 for enumerating the number of code cycles. However, these graphs proved to 
be a useful tool for various problems related to the subject of combinatorics on words 
(e.g. see [3 El It is known that the graph G n can be defined recursively as G n = 
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L(G n -i). The authors were interested in studying other graphs defined recursively using 
the operation of taking line graphs (with natural bases), which could give interesting 
applications. Also, with our choice of graphs (G™ and Kf), it is natural to complete 
them to regular graphs (R% and P") and study these graphs. It turns out that there are 
combinatorial interpretations (relations to other combinatorial objects) for the number 
of independent sets for some of our graphs, and we mention these relations in Sections 01 
and E3 Moreover, we construct a direct bijection describing such a relation for P% (see 
Proposition 13.51 and the discussion that follows). 

2. The numbers g${n). 
Let us first find an explicit formula for g 3 (n). 

It is clear that for any independent set of the graph G%, we can label a node of G% 1 
if this node is in the independent set, and label it otherwise. Thus, our purpose is to 
count the number of triangles having either or 1 in each node and such that no two 
adjacent nodes are both assigned Is. We call such triangles legal. 

In order to get a recursion for g 3 (n), we introduce three auxiliary parameters a n , b n , and 
c n , which are the numbers of legal triangles that, up to rotation, have specific numbers 
in the nodes of the biggest triangle (see Figure EJ). Since we consider only legal graphs, 
the Is in the nodes of the biggest triangle induces Os in certain nodes of a smaller triangle 
(this Os are shown in Figure 0J) . 




o ? o l o i i o i 



Figure 4. Auxiliary parameters a n , b n , and c. 



Considering all the possibilities for the numbers of the biggest triangle, we have that 

03 (n) = #3(n - 1) + 3a n + 3b n + c n , 

where g 3 (n — l) corresponds to all Os, and we have multiple 3 two times because of possible 
rotations. Similarly, we get that a n = g 3 (n — 2) + a n _i, and b = g 3 (n -2). This 

leads to the recursion 

(1) g 3 {n) = 2g 3 (n - 1) + 6g 3 (n - 2) - Ag 3 {n - 3), 
which, under the same initial conditions, is equivalent to the recursion 

(2) g 3 {n)=Ag 3 {n-l)-2g 3 {n-2). 
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We define #3(0) = 1, since we associate the graph G® with the empty graph, in which case 
there is only one independent set, the empty set. Thus, 

g 3 (n) = JL= ((2 + v^ 1 - (2 - v^ 1 ) , 

and the generating function for the numbers gz(n) is 1/(1 — Ax + 2x 2 ). The initial values 
for the numbers gz{ri) are: 

1, 4, 14, 48, 164, 560, 1912, 6528, 22288, 76096, .... 

Preceded by 0, the sequence {g${ri)} is the binomial transform of the Pell numbers 

_ (1 + V2) n - (1 - ygg 

(see [TUl A007070]). These numbers can also be interpreted as maximum bets in a poker 
game (also see [TUl A007070]), where the first player bets 1 dollar into a pot and the ith 
player bets the amount of the (i — l)st player's bet plus the resulting amount of money 
in the pot. Then the number of dollars d n in the pot after n bets is given by 

d n = 2(d n _i + (dn-i - dn-2)) = ^d n ^x - 2d n _2, d = 1, d x = 4, 

which yields d n = g^{n). 

We remark that it would be interesting to obtain recurrence (J2J) directly from the graph 
G3, rather than via recurrence (0). Unfortunately, we were unable to do this. 



3. An algorithm for calculating gi(n), ri(n), kg(n) and pi(n) 

In this section we present an algorithm for calculating gi(n), r^(n), fcg(n) and pe{n) by 
using the transfer matrix method (see [UJ Theorem 4.7.2]). 

3.1. An algorithm for calculating ge(n). In this section, we use the transfer matrix 
method to obtain an information about the sequence of ge(n). 

Similarly to Section [21 for any independent set of the graph G™, consider a labeling 
of G™, where the nodes of of the independent set are labeled 1 and the remaining nodes 
are labeled 0. For a given graph G™, we define the n-th level of G" to be G" \ G™~~ , 
which is isomorphic to G\. Thus, we may think of an independent set of the graph G™ 
as assembled from elements chosen on each level, making sure that when we add a new 
level, we create no conflict with the previous level. 

The collection Ct of possible level labelings is the set of all (0, 1) ^-vectors v = (v±, . . . , ve). 
It will be convenient to define vg + i :— v%. Then v = (v±, . . . , vi) and w = (wi, . . . , we) in 
Ce are a possible consecutive pair of levels in an independent set of G\ (with w following 
v) if and only if 



(3) 



Vi = Vi + \ = or Wi = 0, where i = 1, 2, . . . , £. 
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Thus, to obtain any independent set in the graph G™, we begin with a vector of then 
keep adjoining each next vector w e £i so that it satisfies together with the previously 
chosen vector v G £e, until n vectors have been selected. 

We define a matrix G = Gg, the transfer matrix of the problem, as follows. G is a 
2 e x 2 £ matrix of Os and Is whose rows and columns are indexed by vectors of Li. The 
entry of G in position (v,w) is 1 if the ordered pair of vectors (v,w) satisfies (jHJ), and 
is otherwise. G depends only on £, not on n. Hence, the number of independent sets 
of Gg, ge(n), is the first entry of the vector G n ■ (ui, . . . ,u 2 e) T , where Ui — 1 if there are 
the ith vector v in the collection £g has no two consecutive Is, even after wrapping, (i.e. 
Vi + v i+ i < 1 for all i = 1, 2, . . . , £), and Ui = otherwise. Hence, 

gg(n) = (1,0,... ,0)-G n -( Ul ,...,u 2 ef. 

For instance, when i = 3, the possible level vectors are 

(0, 0, 0), (0, 0, 1), (0, 1, 0), (0, 1, 1), (1, 0, 0), (1, 0, 1), (1, 1, 0), (1, 1, 1), 

except for the last level, where we only have 

(0,0,0), (0,0,1), (0,1,0), (1,0,0). 

If we index the rows and the columns of the transfer matrix G in this order, then we get 



G = 
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The vector (ui, . . . , u 2 i) T in this case is (1, 1, 1, 0, 1, 0, 0, 0). If we now find the first 
entry of the vector (/ — xG)~ x ■ (1, 1, 1, 0, 1, 0, 0, 0) T , where / is the unit matrix, then we 
get that the generating function for g$[n) is given by 1/(1 — Ax + 2x 2 ). We obtain the 
results for larger I similarly. 

Theorem 3.1. The generating functions for the numbers g^n) , g^(n) andg§{n) are given, 
respectively, by 

1 + Ax — x 2 — 2x 3 
1 - 3x - lAx 2 + 15x 3 + 7x 4 ' 
(1 +x)(l + 5x- 8x 2 ) 
1 - 5x - 30x 2 + 69x 3 + 31x 4 - 22x 5 ' 
1 + 10a: - 12x 2 - 50x 3 + 10x 4 + 20x 5 - 12a; 6 
1 - 8x - 66x 2 + 280x 3 + 178x 4 - 532x 5 - 84x 6 + 108x 7 ' 

We remark that the algorithm for finding the generating function for gg(n) has been 
implemented in Maple, and yielded explicit results for £ < 6. 
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3.2. An algorithm for calculating r^(n). In this section we use the transfer matrix 
method to obtain an information about the numbers rp(n). This case is similar to that of 
gt{n) with some small differences. 

We partition i?" into levels just as in the case of G\, so the ra-th level of i?" is . 

The collection of possible levels Li is the set of all ^-vectors v of Os and Is such that 
there no consecutive Is in v, that is, t>j + v i+ i ^ 2 (where we define vp + i := Vi). Clearly, 
the set Ce contains exactly Lg vectors where Li is the £th Lucas number. For instance, 
C 3 contains the vectors (0, 0, 0), (0, 0, 1), (0, 1, 0), and (1, 0, 0). 

The condition that vectors v and w in Ce are a possible consecutive pair of levels in an 
independent set of i?" is given by (jHJ) just as for G™. To obtain any independent set in 
the graph i?™, we begin with a vector of Ci, then keep adjoining each next vector w G £e 
so that it satisfies Q together with the previously chosen vector v G Ci, until n vectors 
have been selected. 

We define a matrix the transfer matrix of the problem R = Ri in the same way as 
G = Gi in subsection 13.11 Then R is an Lg x Lg matrix, and the number of independent 
sets of R™, ri(n), is the first entry of of the vector R n ■ 1 where 1 = (1, 1, ... , 1). Hence, 

r,(n) = (l,0,...,0)-iT-l. 

For instance, when I = 3, the possible level vectors in an independent set are 

(0,0,0), (0,0,1), (0,1,0), (1,0,0). 

If we indexed the rows and columns in this order, then the transfer matrix is 

"l 1 1 l" 

10 1 
110 0' 
10 10 

If we find the first entry of the vector (/ — xG)^ 1 ■ 1, we get that the generating function 
for r^{n) is given by 1 2 . The initial values for the numbers r 3 (n) are 

1, 4, 10, 28, 76, 208, 568, 1552, 4240, . . . 

This sequence appears as A026150 in [TOj . 

Similarly to the case i = 3, we obtain the following results for £ = 4,5, 6. 

Theorem 3.2. The generating functions for the numbers r^n) , r 5 (n) andr§{n) are given, 
respectively, by 

1 + Ax - Ax 2 
1 — 3x — Ax 2 + Ax 3 ' 

1 + 7x - Qx 2 
1 — Ax — 8x 2 + 6a; 3 ' 
1 + 12x - 2Ax 2 + 8x A 
(1 - 8x + Ax 2 + Ax 3 ) (1 + 2x - 2x 2 ) ' 
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We remark that the algorithm for finding the generating function for r^(n) has been 
implemented in Maple, and yielded explicit results for i < 6. 

3.3. An algorithm for calculating k^{n). In this section we use the transfer matrix 
method yet again to obtain information about the sequences k^{n). This case is also 
similar to that of ge{n), so we will only sketch it briefly. 

We partition into levels just as in the case of G™, so the n-th level of is K™\K™~ X . 
The collection of possible levels Li is the set of all ^-vectors v of Os and Is. Vectors v 
and w in Li are a possible consecutive pair of levels in an independent set of K™ if they 
satisfy ©• 

We define the transfer matrix of the problem, K = K#, in the same way as G(_. Then K 
is a 2 e x 2 e matrix, and the number of independent sets of K^, ki(n), is the first entry of 
of the vector K n ■ (ui, . . . , u 2 e) T where «j = 1 if the zth vector in the collection Ce contains 
at most one nonzero entry. Hence, 

k e (n) = (l,0,...,0)-K n -( Ul ,...,u 2e f. 

For instance, when i = 3, the possible level vectors in an independent set are 

(0, 0, 0), (0, 0, 1), (0, 1, 0), (0, 1, 1), (1, 0, 0), (1, 0, 1), (1, 1, 0), (1, 1, 1). 

If we index the rows and columns in this order, then the transfer matrix is 



K = 
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If we find the first entry of the vector (J — xK)~ l ■ (ui, . . . ,u 2 e) T , we get that the 
generating function for k^{n) is given by 1 _ 4a; 1 +2x2 . In particular, we get that fc 3 (n) = g^{n) 
which is can also be seen directly from the definitions. Similarly, we have the following 
result. 

Theorem 3.3. The generating functions for the numbers k^{n), k 5 (n) and k 6 (n) are 
given, respectively, by 

1 + 2x + 3x 2 
1 - 3a; - Ux 2 + 15x 3 + 7x 4 ' 
1 + x + 12x 2 - 8x 3 
1 - 5x - 30a; 2 + 69a; 3 + 31x 4 - 22a; 5 ' 
1 - x + 38a; 2 - 72a; 3 - 8x 4 + 30a; 5 
1 - 8a; - 66a; 2 + 280a; 3 + 178a; 4 - 532a; 5 - 84a; 6 + 108a; 7 ' 

We remark that the algorithm for finding the generating function for k^{n) has been 
implemented in Maple, and yielded explicit results for £ < 6. 
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3.4. An algorithm for calculating pe(n). In this section we use the transfer matrix 
method to obtain information about the sequences pe(n). 

In this section we use the transfer matrix method once more to obtain information 
about the sequences fc^(n). This case is also similar to that of ge(n), so we will only sketch 
it briefly. 

We partition P" into levels the same way as G™, so the n-th level of P™ is P"\P™ -1 . 
The collection of possible levels Ce is the set of all ^-vectors v of Os and Is. Vectors v 
and w in Ce are a possible consecutive pair of levels in an independent set of P™ if they 
satisfy (J3J). 

The collection of possible levels Ce is the set of all £- vectors v = • • • , Vi) of Os and 
Is such that v\ + • • • + vi < 1. Clearly, the set Ce contains exactly £ + 1 vectors which 
are (0, • ■ • ,0) and (0, . . . , 0, 1, 0, . . . , 0). The condition that vectors v and w in Lt are a 
possible consecutive pair of levels in an independent set of P" is given by (JHJ). 

We define the transfer matrix of the problem, P = Pi, in the same way as Ge, Re and 
Kp. We define a matrix P = Pi, the transfer matrix of the problem, as follows. P is an 
{£ + 1) x (£ + 1) matrix of 0s and Is whose rows and columns indexed by vectors of Ci. 
Therefore, it is easy to see that 



1 1 




1 
1 





1111 
111 



1 
1 1 



(e+i)x(e+i) 

Theorem 3.4. The generating function for pe(n) is given by 

l + 2x 



£ 

n>0 



p e {n)x 



1 - (£ - l)x - 2x 2 ' 



Proof. We want to find the first entry of the vector (J — rcP) -1 ■ 1 which means we must 
find the first row, say (ex, ... , e^+i), of the matrix (J — xP)^ 1 . By solving the system of 
equations 



(/ - xP) 



-i 



ei, 



e m ) T = (l,0,...,0) T 



we get that 



ei 



I -{I- 2)x 



and ej 



x 



l-(£-l)x-2x 2 J l-(£-l)x-2x 2 
Hence, the first entry of the vector (J — xP)^ 1 ■ 1 is given by 

i+i 



for j > 2. 



l + 2x 



8=1 



- 1 )x-2x 2 



□ 
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For instance, when t = 3, the generating function for Pz{n) is given by ^-vx-vx 1 • ^ ne 
can see, in particular, that pzin) = rs(n), which follows directly from the definitions. 

In the case i = 4 the initial values of the numbers p±{n) are 

1, 5, 17, 61, 217, 773, 2753, 9805, 34921, 124373, . . . 

The same sequence turns out to appear in jH] (see [HI A007483]). Thus, the following 
proposition is true: 

Proposition 3.5. The number of independent sets in the graph P£ is equal to the number 
of (possibly empty) subsequences of the sequence {1,2, ...,2n + 1} in which each odd 
member has an even neighbor. 

Here, the neighbors of an integer m are m— 1 and m+1. In the case n = 1, the sequences 
appearing in the proposition are e, 2, 23, 12, 123, where e is the empty sequence. In this 
case, we can find a direct bijection between the objects in Proposition 13.51 as described 
below. 

We start by labeling the vertices of the (innermost) level n of P% clockwise by 2, 23, 12, 
123. The level n — 1 is labeled starting from the vertex immediately to the left of vertex 
labeled 2 as follows: 4,45, [3]4, [3]45 (the meaning of brackets will be discussed below). 
More generally, for i < n, given a level n — i + 1 labeled clockwise with 2i, 2i{2i + 1), [2i — 
l]2i, [2i — l]2i(2i + 1), we label the (next outer) level n — i clockwise from the inside out 
with 2i + 2, (2i + 2) (2i + 3) , [2i + 1] (2i + 2) , [2i + 1] (2i + 2) (2i + 3) starting from the vertex 
immediately to the left of vertex labeled 2i. (See Figure El for the case n = 3.) Each 
independent set has at most one vertex on each level. Now, given any independent set in 
P£, we can write the labels of its nodes in increasing order and delete any integer [2i — 1] 
in brackets if the sequence also contains 2i — 2 or 2i — 1 without brackets. Erasing all 
brackets now, if any, we obtain a sequence from Proposition 13.51 




Example 3.6. The independent sets {45, [5] 67}, {4, [5] 67}, {12, [5]6}, {[3]4, 67} correspond 
to the sequences 4567,467,1256,3467, respectively. 
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For convenience, we will write down the set of nonadjacent labels at level n — i for each 
label at level n — i + 1. 

(no label) e ■-> (2i + 2), (2i + 2)(2i + 3), [2i + l](2i + 2), [2i + l](2i + 2)(2i + 3) 
2i i — * [2» + l](2i + 2), [2i + l](2i + 2)(2i + 3) 
(4) 2i(2i + 1) ^ (2i + 2), [2i + l](2i + 2)(2i + 3) 

[2i - 1]2» i — ► (2i + 2), (2i + 2)(2i + 3) 
[2i - l]2i(2i + 1) i-> (2i + 2)(2i + 3), [2i + l](2i + 2) 

Now it is not difficult to construct an independent set given a sequence of Proposi- 
tion 13.51 We partition the sequence of integers from 1 to 2n + 1 as follows: 

123 | 45 | 67 | ... | 2n(2n + 1), 

then choose the vertices of the independent set in the order of increasing labels using the 
rules Notice that the label of the vertex at level n — i + 1 must contain 2i. 

Example 3.7. 

4567 ^ 123|45|67 i-> ( e ,45, [5]67) i-> {45, [5]67} 
467 i-> 123|45|67 ^ (e, 4, [5]67) i-> {45, [5]67} 
1256 i-> 123|45|67 ^ (12, e, [5]6) i-> {12, [5]6} 

It can be shown that the two maps described above are inverses of each other based on 
the recursive structure of sequences under consideration. 
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